Probing the CP-even Higgs Sector via H 3 -> H 2 Hi in the Natural NMSSM 



o 

(N 
4=> 

IX, 
(N 



43 

Oh. 
i , 

Ph. 

<D 

43 



> 

m 
in 

q 
o 



x 

S3 



Zhaofcng Kang, 1 '* Jinmian Li, 2 ' ^ Tianjun Li, 2 ' 3,4 '- 1 - Da Liu, 2,§ and Jing Shu 2 ' 11 

1 Center for High- Energy Physics, Peking University, Beijing, 100871, P. R. China 
State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 
Chinese Academy of Sciences, Beijing 100190, P. R. China 
'School of Physical Electronics, University of Electronic Science and Technology of China, Chengdu 610054, P- R- China 
* George P. and Cynthia W. Mitchell Institute for Fundamental Physics and Astronomy, 
Texas A&M University, College Station, TX 77843, USA 
(Dated: February 25, 2013) 

After the discovery of a Standard Model (SM) like Higgs boson, naturalness strongly favors the 
next to the Minimal Supersymmetric SM (NMSSM). In this letter, we point out that the most 
natural NMSSM predicts the following CP-even Higgs Hi sector: (A) H2 is the SM-like Higgs boson 
with mass pushed-upward by a lighter Hi with mass overwhelmingly within [mH 2 /2,mH 2 ]; (B) 
ms 3 — 2/it/ sin 2/3 > 300 GeV; (C) H3 has a significant coupling to top quarks and can decay to 
H1H2 with a large branching ratio. Using jet substructure we show that all the three Higgs bosons 
can be discovered via gg — ► Hz — » H\Hi — > bblujj at the 14 TeV LHC. Especially, the LEP-LHC 
scenario with Hi ~ 98 GeV has a very good discovery potential. 

PACS numbers: 12.60.Jv, 14.70.Pw, 95.35.+d 



Introduction: Supersymmetry provides the most ele- 
gant solution to the gauge hierarchy problem in the Stan- 
dard Model (SM). In the supersymmetric SMs (SSMs) 
with i?-parity, we can not only achieve the gauge cou- 
pling unification, but also have a cold dark matter candi- 
date. Recently, the discovery of a SM like Higgs boson at 
the LHC with mass mh around 126 GeV [1] has deep im- 
plications to the SSMs. Although such a relatively heavy 
Higgs boson mass can be achieved in the Minimal SSM 
(MSSM), it generically incurs a large fine-tuning (For the 
possible solutions, see [2]). By constrast, the next-to-the 
MSSM (NMSSM) with an extra SM singlet Higgs field 
S is strongly favored by naturalness [3], due to origi- 
nally its dynamically solution to the Higgs bilinear mass 
/i problem and now the SM-like Higgs boson mass en- 
hancement via the relatively large Higgs trilinear Yukawa 
coupling A in the superpotential and singlet-doublet mix- 
ing effect [4-7] . The natural NMSSM may leave hints at 
the light stop sector, but the search is rather model de- 
pendent [8, 9] and barely has relation with Higgs sector 
(Recent attempt to search for the light stop utilizing the 
properties of the SM-like Higgs boson was done in [10].). 

In the natural NMSSM, the second lightest CP-even 
Higgs boson H2 is indentified as the SM like Higgs boson, 
while the lightest CP-even Higgs boson Hi has dominant 
singlet component. Thus, the H2 mass can be pushed- 
upward via the singlet-doublet mixing effect [4-7] . Such 
a scenario can explain the possible di-photon excess from 
Higgs decays [4, 11, 12] since the significant mixing effect 
reduces the decay width of H% —> bb and the light charged 
Higgsino may increase the Higgs decays to diphotons. In- 
terestingly, Hi may be used to interpret the slight LEP 
excess for the Higgs mass around 98 GeV [14] (It receives 
some interest [15, 16] recently.), or the LHC excess for the 
Higgs mass around ~113 GeV [17]. A scenario with two 



light higgs and a low-mass pseudoscalar in NMSSM has 
been discussed in [18]. More noticeable features emerge 
when we take the heavy CP-even Higgs boson H3 into 
account. In this letter, we consider the CP-even Higgs 
sector in the natural NMSSM. We point out that natu- 
ralness implies the H3 mass range m# 3 G [300, 600] GeV 
and its significant triple Higgs coupling with Hi and Hi- 
Such a Higgs sector structure leads us to investigate the 
discovery potential of the whole CP-even Higgs bosons 
from the process gg — > H3 — > H1H2. With jet substruc- 
ture, we show that all three CP-even Higgs bosons Hi 
can be probed at the 14 TeV LHC. Our search strategy 
is specially suitable for the LEP-LHC Higgs bosons but 
also applies to the general pushing-upward scenario. 

Light Higgs Bosons in the Pushing-Upward Sce- 
nario: The SM-like Higgs boson can be accommodated 
without recurring severe fine-tuning, and we can show 
that the whole Higgs sector is light. Restricted to the 
Z3— NMSSM, naturalness conditions point to a predic- 
tive parameter space 



A : 0.6-0.7, tan/3 : 1.3-3.0, 
H = Xv s : 100 GeV - 200 GeV, 



(1) 



where tan (3 is the ratio of the vacuum expectation values 
for two Higgs doublets, and K is the singlet cubic coupling 
in the superpotential. Also, k is constrained by perturba- 
tivity, and typically is no more than half of A. The stop 
sector should be sufficiently light, e.g., mj L = mj R = 500 
GeV, and a flavor safe choice A t = —500 GeV. Their 
concrete values will not qualitatively affect our following 
discussions. 

Importantly, A\ can be further determined in the 
pushing-upward mixing scenario. The Higgs mass square 
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matrix in the Goldstone basis is 



(M|)n = M 2 A + (m| 
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where Ml 



2Av s (A\ + kv s )I sin 2/3 defines the largest 
scale among these elements. Let the orthogonal matrix 
diagonalizing Mj be O: T D\&g(m? H ,fn 2 H2 ,m% 1 )0 = 
Mg. The singlet-doublet mixing effect can be approx- 
imately studied by decoupling the entries involving the 
first state. Rcf. [4] found that, in the case with a large A 
and small fj,, the realization of pushing-upward scenario, 
which requires (Mg)^ < (AfJ)22, necessitates a cancela- 
tion to reduce the large non-diagonal element (MD23: 

1 — (Ax/2/j, + k/X) sin 2/3 ~ 0. (3) 

Thus, A\ is largely determined by fj, and tan /3, and to a 
less degree, by k. Then we have 
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2/< 



K sin 2/3 
~ A 2 



(4) 



, sin 2/3 , 

Recall that k < A, so, to a good approximation, we get 
rriH 3 ~ Ma — 2fi/ sin 2/3, which is about 2.5//, relating 
the H3 mass directly with the weak scale naturalness. 

We now summarize the Higgs spectra in the natural 
NMSSM under consideration. First, all the Higgs fields 
are properly light. H3 and its SU(2)l partners, the 
charged Higgs bosons i? ± and the heavy CP-odd Higgs 
A2, take roughly degenerate masses Ma- Hi is SM-likc 
while Hi is even lighter. Hi is a SM singlet like and then 
can be allowed by the LEP experiment. Note that mH x 
is most likely to fall into the region [m^/2,m/j] with the 
lower bound set by forbidding the decay Hi — > H1H1 
(Ref. [19] considered such case.). Otherwise it tends to 
be the dominant decay mode of H^- In addition, the 
lightest CP-odd Higgs boson Ai also has a mass around 
the weak scale. Moreover, a pair of charginos and three 
neutralinos, consisting of the Higgsinos and singlino, are 
light as well. All of them may be detectable at the LHC 
and here we focus on the CP-even Higgs bosons. 
Hi— couplings: The Higgs signals at colliders are sensi- 
tive to their mixing angles whose effects, in a standard 
form, are described by the tree-level Lagrangian: 



Dn , z ^H l ZZ + r t , w ^k Hl W + W- 
V2v v 



V2i 



HijkHiHjHk, 



(5) 



with v pa 174 GcV. r^y. etc., encode the deviations of 
Hi from hsM- For instance, we have 



ri.v = O32, r 2 ,v = O22, T3,v = Oi2 



(6) 



We also include the triple Higgs couplings, which will 
play a crucial role in the search for Higgs bosons. 

We now present the features of H 3 couplings. Firstly, 
note that (M|) 12 is a small entry and we can express it in 
terms of O and rn 2 H . . Since m# 3 is a few times of m^f 2 1 , 
then it is not difficult to obtain the upper bound 



Oi 



so, < (Mj) 12 /m^ 3 ~ {Mi,)i 2 /(Mi,)i 



(7) 



where (M|)n gives the dominant contribution to m^. 
Therefore, the trilinear couplings between H 3 and the 
weak gauge bosons are negligibly small. Next, the re- 
duced couplings of H3 to the bottom and top quarks are 
given by 



C.3,5 = - On tan/3 + 0i2 w -On tan/3, 
C3 1 =Oii cot /3 + 0i2 w On cot /3. 



(8) 



Owing to a relatively small tan /3 in the natural NMSSM, 
H3 coupling to the bottom quark is not enhanced while 
its coupling to the top quark is significant. They have 
crucial implications to the collider phenomenology of H3, 
e.g., it can be considerably produced at the LHC by 
virtue of the significant coupling to gluons: 



3 ,s =l.O30 2 ,t - O.O60 2 , fc « On cot /3. 



(9) 



Finally, the triple Higgs coupling H3H2H1 receives two 
possible large contributions and is given by 



^123 



A^a 
V2 
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A A\ 



XA X 
V2 



(10) 



It thus has a large \A\ enhancement and leads to H3 — > 
H1H1 decay width at the GeV scale and dominates the 
H3 Higgs-to-Higgs decay, as provides the most promising 
discovery prospect for H 3 and Hi, similarly to Ref. [13]. 

We now turn our attention to the lightest Higgs bo- 
son Hi. Interestingly, the LEP collaboration reported 
(with an signal significance 2.3 a) a slight excess of events 
for a Higgs boson with mass ~ 95 — 100 GeV [14]. Al- 
though our discussions on the Higgs bosons and the en- 
suing search strategy are not restricted to this case, it is 
tempting to interpret Hi as the source of this excess. So 
we have 



r 2 



Br(i?i -> 66) 
Br SM (iii ->• bb) 



0.1-0.25. 



(11) 



For ma S 100 GeV, its decay to 66 nearly determines 
its total width. Thus, the LEP requires C^yv ~ 0.3 
which is a typical value expected from the mixing Higgs 
sector. 
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FIG. 1: A plot on the aH 3 -m,H 3 plane, with color code de- 
noting rriHx ■ Large inverted triangle points satisfy the LEP- 
LHC scenario. We use NMSSMTools 3.2.1 [20], set ^ C 
[300, 500], ^ C [-300, 0], and require ^ C [125,127] 
and signal strengthes R2, 99(77) C [1.4,1.6], R2, gg {VV) C 
[1.0,1.3]. 



Signature and backgrounds: In light of the previ- 
ous analysis, the signature gg — > H3 — > H\{— > bb)H 2 (-^ 
WhWt) is promising, where we denote Wh as hadronic 
decaying W boson and denote We as leptonic decaying 
W boson. The We will suppress the enormous QCD back- 
grounds. The total cross section is 



=0.2 



Br(H 3 -> H1H2) Br(#i -> bb) 



0.4 J ' 20% 
Br{H 2 -> W e W h ) a GF (hs M ) 



90% 



28% 



10 pb 



pb, 



(12) 



where £ = e, ji. The numerical results are shown in Fig. 1, 
where a plot of the distribution of <jh 3 on the m# 2 — 
plane is presented. It can be seen that its values cluster 
well for a given m# 3 (typically within only a few times) , 
in particular for heavier H3. 

We implement the simplified model for Higgs bosons in 
Feynrules [21] to generate the UFO format of the effec- 
tive model for MadGraph5 [22], where the parton- level 
signatures are generated. 

The semi-leptonic ti pair production is the dominant 
background (BG), with the NNLO cross section rj 240 
pb [23]. The subdominant BG We + 66+jets has cross 
section depending on the renormalization scale, roughly, 
about 40 pb. Other backgrounds can be neglected in 
our signal region. BGs are generated using MadGraph5. 
To avoid double counting, we adopt the modified version 
of MLM-matching [24] with xqcut = 15 GeV. For the 
latter BG, we include up to 2 additional jets and set the 
fc— factor to be 2. 

We use PYTHIA6.420 [25] for decaying particles, 
parton-showering and hadronization. However, in or- 
der to employ the BDRS procedure later, we turn off 
the B— hadron decays in Pythia. The produced objects 
are then converted to the HepMC [26] event format and 
passed to Fastjct 3.0 [27] to cluster the final states. The 
final visible particles arc requited to have pr > 0.1 GeV 



and 1 77 1 < 5.0 which are defined as tracks hereafter. Lep- 
tons from signal events should be isolated, otherwise they 
are combined with the tracks to reconstruct fat jets later. 
Additionally, signal leptons arc required to have < 2.5 
and pt > 10 GeV. We take b— tagging efficiency of 70% 
with the other light quark mis-tagged probability 1%. 

We choose the C/A algorithm [29] with radius R=1A 
and pt>40 GeV to cluster the tracks and form fat jets. 
Following BDRS [28], we first break the hard fat jets into 
subjets ji t 2 with masses rrij 1 2 . Next, a significant mass 
drop rrij 1 < fi rrij with [1=0.667 and not too asymmet- 
ric splitting, i.e., y =min(p|, Ji , p^ h )AR'j uj2 /m 2 J > y cut 
with j/ cut =0.09 (Ai?| - 2 is the angular distance), are re- 
quired. If the above criterion are not satisfied, we will set 
j = ji and go back to decomposition. Finally, we filter 
the Higgs neighbourhood, resolving the fat jets on a finer 
angular scale Rfnt = min(0.35, Rj 1 j 2 /2) and taking the 
three hardest objects, with the remains identified as the 
underlying events contamination and hence dropped. 
Events selection and results: Two basic cuts are im- 
posed to trigger our events. Firstly, at least two filtered 
fat jets are required. One of them has two leading subjets 
which pass &— tagging and satisfy \rj\ < 2.5, and then is 
identified as the H±— jet. Among the remaining fat jets, 
the one with highest pt is regarded as the Wh— jet [30]. 
Secondly, the events must contain exactly one isolated 
lepton. 





FIG. 2: Distribution for trigged signal and background of 
Pt,h x , m Hl , rnjjw, Afaj, AR Hlb i, Mc- The number of events 
have been normalised to 14 TeV 500 fb^ 1 , and the signal is 
400 times amplified. 



For illumination, here we will take a benchmark point 
inspired by the LEP-LHC Higgs scenario: mjjj = 98GeV, 
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m H2 = 125 GcV as well as m H3 = 400 GeV. Figure 2 
shows the distributions of some important kinematic vari- 
ables. In terms of the plots, we display the cut flow: 

• Cutl: The relatively large mass splitting between 
H3 and Hi gives H\ a boost. Therefore, we require 
p Tbb > 150 GcV, prjjlv > 120 GcV, and \p T bb - 
p T ,jj£u\ < 20 GcV. 

• Cut2: It is observed that the longitudinal momen- 
tum of the neutrino from W decay is gcncrically 
small, and hence mjy 2 3 can be approximately re- 
constructed by assuming p z>v — 0. Practically, cuts 
based on this assumption are sufficiently good. So 
we impose: 95 GeV < mjjj < 100 GcV, rrijjjh, < 
150 GeV, and m bb]]lv < 440 GeV. 

• Cut3: Because H2 has spin-0 and W only couples 
to the left-handed fermions, the lepton from We 
will align with one of the jets from Wh decay. It 
allows us to impose a cut |A0^ | < 1.5, namely 
the azimuthal angles difference between the signal 
lepton and (one) jet being sufficiently small. 

• Cut4: The filtered -Hi-jet actually contains three 
subjets, the bb and a radiated gluon. So the -Hi-jet 
and its bb subsystem must have a very small angle 
distance. By contrast, the angle distance between 
the .ffi-jet and WVjet is much larger. Thus, we 
require AR Hi bb < 0.01, and 2.6 < AR HlWh < 3.4. 

• Cut5: We also impose the cluster transverse 
mass of decay product of the Hi'. Mc = 

\jPTjH + m %i +$t< 220 GcV. 

With the above cuts, we obtain the signal significance 
4.42 a excess for the LEP-LHC benchmark point at 14 
TeV 500 fb~ x . The cut efficiency and the signals are 
presented in Table I. 





ti 


W(- 


-» lvjj)bb + jets 


Signal 


Total 


1.2 x 10 8 




1.91 x 10 7 


1.25 x 10 4 


Trigged 


4.95 x 10 6 




1.45 x 10 6 


1456.75 


Cutl 


3.77 x 10 5 




1.61 x 10 5 


639.5 


Cut2 


1932 




203 


119.75 


Cut3 


1512 




155.2 


105.5 


Cut4 


108 




47.75 


56.25 


Cut5 


84 




47.75 


55 



TABLE I: Number of events after each cut for background 
and signal (normalized to 500 fb^ 1 ). The signal significance 
SyV S + B has reached to 4.02 and with the precise 4.42 a 
excess for the LEP-LHC benchmark point. 

Since Fig. 1 shows obvious cluster behavior, the whole 
parameter space with pushing-upward effect can be ex- 
plored. Using BDT analysis [31], we consider six repre- 
sentative points to demonstrate the search prospect, and 



the discovery signal significance for each case is given 
in Table II. Some observations can be made: (A) For a 
given , a lighter Hi shows better discovery potential; 
(B) Increasing H 3 mass helps to boost Hi but the cross 
section is reduced. Thus, a moderately heavy H 3 ~ 400 
GeV and relatively light Hi have the most promising dis- 
covery potential; (C) Most of the parameter space is dis- 
coverable except for simultaneously light H3 and heavy 
Hi, e.g., the benchmark point Bl, despite of a rather 
large cross section, has a quite low signal significance. 

The situation can be further improved when we take 
i? ± into account. _ff ± can be produced associated with 
a single top, with a moderately large cross section at the 
small tan/3 region. Moreover, it can decay to Hi and 
W with a substantial branching ratio and hence provide 
a way to probe Hi and i? ± . In this case H^ can be 
produced with a larger px and the boost will be easier for 
lighter H , not very sensitive to m^. So it can provide 
a complementary or even more promising channel for the 
pushing-upward scenario. We leave it for the future work. 





m Hl (GeV) 


m H3 (GeV) 


a (fb) 


s 




VS+B 


Bl 


100 


300 


70 


0.81 


B2 


65 


300 


50 


3.84 


B3 


98 


400 


25 


4.73 


B4 


65 


400 


20 


7.68 


B5 


100 


600 


2 


2.79 


B6 


65 


600 


2 


4.99 



TABLE II: Discovery signal significances for 6 representative 
points at 14 TeV 500 fb^ 1 . We design 25 kinematic variables 
for BDT analysis [31]: $ T , Pt,w, m w , n jet , Pr,b 1 , Pr,b 2 , Pr,e, 
mT,iv, Pr,wji, Pt,wj 2 , PT.jjtv, ARij, A(j>i w , pt,h 3 , m-iv, E\ v , 
and m,H 3 ■ 

Conclusion: We pointed out the specific features in the 
CP-even Higgs sector of the natural NMSSM, and showed 
that all three CP-even Higgs boson Hi can be probed at 
the 14 TeV LHC 
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